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THEOl^TICAL ANALYSIS OF AN AUGMENTOR WING 
FOR A VTOL FIGU'IER 


by Mnrnix P. B, Dillcnius and Michaol R. Mendenhall 

SUMMARY 

An analytical method has heon developed foi* predicting forces and 
moments acting on augmontor wings in hover or forward flight for pre- 
scribed ogcotor jot characteristics. I’ho method is based on incorapros- 
siblo, potential flow theory. Attached flow is assumed, Two potential 
flow models are incorporated in the prediction metliod: a model for the 

wing/flap system and a model for the jet and its v/ako. A tlu-oo-dimen” 
sional nonplanar vorto.x lattice is laid out on the chordal planes of the 
augmontor wing components. Flap surfaces and jot can bo at largo angles 
to the oncoming stream, ilorsoshoo vortex strengths are determined from 
th<2 application of the flow tangency condition at points on the surfaces 
including intorfaranca effects from tlic jet. The jet is made to expand 
from the primary no' 45 ?les to the diffusor exit, and a distribution of vor~ 
ticity is placed on the jot boundary to model entrainment, in addition, 
blockage panels are positioned on the boundary of the jot wake downstream 
of the augmontor wing. Jot wake contorline location and spreading rate 
are taken from exporimontal data. The solution proceeds in an iterative 
manner using the two flow models in scguonco and comparing the predicted 
diffusor exit velocity with the specified velocity. 

Experimental data involving VTOL fighter-type augmontor wings are 
not readily available; however, some comparisons are shown using data on 
a V/STOL transport model in forward flight. Good agreement is shown at 
low thrust settings. For high thrust settings, the present method over- 
predicts tl20 overall lift and drag and undorprodicts the pitching moment. 
The discrepancy is at least partly attributable to the attached flow 
assumption. 



INTHODUCl’ION 


1‘hoKa has boon on .ineroasing inUoKost in tho clevolopmont of VTOL 
high porfomanco aircraft;, Ono Gonaopb usotl to achiovo VTOL capability 
is tho nugmontor wing. In this concopt» tho propulsion system in into- 
grated with the lifting surfneos for generating direct lift. The aug- 
roontor wing of interest hero consists of a fixed wing with a sot of 
largo flap surfaces at tho trailing edge which can bo deflected at large 
nngloa to guide and control the engine exhaust gas. The gas or primary 
jot emanates from nonzles at tho trailing edge of tho center flap 
positioned above tho forward and aft flaps designed to form a diffusor 
or ejector system, in this process,^ the jot entrains secondar'y air 
drawn into tlio diffusor and mixes with it resulting in an augmentation 
in a\ass flow. As a result ,» the overall thrust is larger than tho primary 
thrust generated by the engine e.xhaust dviverted from tho engine to tho 
augmentor wing. FurthaA’morOi additional lift is generated by the wing 
and flap system due to interferonco offocts induced by tho augmented 
jot. flOV'?. 

The objeetivo of this report is to describe an analytical method 
developed for tho prediction of the external aoi'odynamie elwracteristias 
of an augmentor wing for specified jet velocities at the primary jet 
noarslo and tho diffusor exit. Plow conditions include )iovor and forward 
flight. This effert is a necessary precursor to the prediction of tho 
longitudinal aerodynamic forces and moments acting on a complete VTOL- 
typo fighter aircraft with one or more augmentor wing systems. The work 
performed during this phase was funded jointly by NASA and NAVAIh. 

Several analytical methods based on jot flap theory liavo been 
developed to analyze wind-on porCox-manco of augmentor wings, Two-dimoii’* 
sional approaches are given in rofarencos 1 through 4. The thin jot 
flap theory is not a realistic model for tho "thick" augmentor wing 
jots. Bovilaqua, roference S, presents a two-dimensional analysis for 
tho static (wind-off) case using a viscous solution for the flow through 
the diffusor and on inviscid solution for tho external flow. Tho viscous 
inner flow solution is based on a turbulence kinetic energy model to 
account for mixing. Very recently, the inner solution was extended to 
account for hyper-mixing nozzles (ref. 6). Bevilaqua's methods can pre- 
dict details of the jet flow inside the diffusor including average jot 
velocities at tho exit of the diffusor. 
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The pircsont roothod is direcUod towai‘ds tho throc-diinonsional case 
of a finito-?ispoGt-iratio wing with swoop and taper having a highly 
dofloctod, thick »augmontor jot. It is dosirablo, as well, to account 
for spanwiso variation of tho jot widtJi and jot volocitioa for a 
realistic throQ-diJiionsional easo. The approach Gonsistn of dovoloping 
models for tho wing-flap system and those aspoeta of the augmontoi* jot 
that aro required to properly model tho intoi^foronGo effects of the jot 
on the wing-flap loading . The details of the flow inside the auginontor 
and tho mixing of the primary and secondary flows are not addressed, 
and the augmontor jot eharaetoristics are assumed to be prescribed 
initially in terms of parameters which yield the primary jot velocity 
and tho mixed flow velocity at tho ejector exit. Potential flow methods 
are employed, and no flow separation on tho wing or flap surfaces is 
assumed to occur , 

Recently, a Navy fighter prototype employ ing an augmontor canard 
and wing was built and tested, it is dosignated XFV-12A and is described 
in reference 7. No tost data is readily av^ailablo for this configuration, 
t'lie basic features, including the augmontor wing systems, aro indicated 
in figure 1. In order to evaluate the methods and provide some com- 
parisons witli data, calculations were made on a rectangular, finite- 
aspoct-i‘atio V/STOb transport configuration for whicli wind tunnel data 
are available (rof. 8). Unlike tho fighter configuration, the transport 
wing has a short chord and high thickness ratio and tho flap segments 
also have high thickness ratios . A ci*oss section of the configuration 
is shown in figure Results are presented on this configuration, and 
discussion is presented on the flow and load characteristics as illumi- 
nated by the comparisons between prediction and data. 
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nspQct ratio 

width of jet measured perpendicular to jet center 

axial force coefficient, axial forco/qSj^gp, figure 3 

pitching moment coofficiont, pitching moment/qSppp-Gppp, figure 3 

normal force coefficient, normal forco/qSppp, figure 3 

jet momentum coefficient, cqn. (9) 

average wing chord 

aerodynamic force vector 

roforcncQ length, " c 

mass flow rate 

jet mass flow rate at jet exhaust 

entrained air mass flow rate 

vector normal to surface 

dynamic pressure 

resultant flow velocity vector 

length along jet centerline 

reference area 

thrust force 

free stream velocity 

flow velocity 

wing coordinate system, figure 3 
angle of attack, degrees 
flap deflection angle, degrees, figure 2 
vortex vector with strength V 
vorticity strongth/unit length 
thrust augmentation ratio 

nondimens ional length along jot conterlino 
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SYMBOLS (Concluded) 


contor of 0Yavity 
jot ccntorline 
diffusor exit 

initial (at cantor flap exhaust nozalo location) 



TECHNICAI, Al>PROACU 

The basic objccbivo of the analytical method is to represent the 
solid surfaces of the auejmentor wing by a distribution of lifting-typo 
singularities for the purpose of calculating aerodynamic loads including 
jet intorforonco effects. To this endf the jot must bo represented by 
suitable singularities from its origin to some distance downstream of the 
diffusor in order to model the jet induced effects. A nonplanar vortor 
lattice is laid out on the augmontor wing components, This lifting sur- 
face moecl is basically a modified form of the vortex lattice approach 
described in roforonce 9. The vortex distribution used to model the 
jet over its length is an adaptation of the jot model used in the work 
described in references 9 and 10, In addition, the jet wake (jot dovm- 
stroam of the diffusor exit) makes use of blockage panels on its bound- 
ary. This part of the jot model is an extension of the blockage panel 
scheme described in reference 11, 

The lifting surface and jot models have boon programmed and arc 
used in sequence to form an iterative approach to tho calculation of the 
longitudinal oorodynaraic charocteristics (lift, drag, pitching moment) 
of augmontor wing configurations. In this process, tho assumption is 
made that flow ontrainmemt by tho primary jot is sucli that tl^o flow at 
the diffusor exit is completely mixed. This condition, to a first 
approximation , is approached in practice when hyper-mixing primary 
nozzles are employed in conjunction with limited boundary layer blowing. 

in tho following sections, brief descriptions are given of the lift- 
ing surface and the jot flow models. Attention will be focused on tho 
special features developed to apply the floxi/ models to augmontor wings. 
Tho iteration scheme is described in the section concerned with tho 
method of solution. 

Vortex-Lattice Model for Augmontor Wing/Flaps System 

General description .- The vortex-lattice lifting-surface model used 
to represent the solid surfaces is a modified version of the model used 
in reference 9. Configurations of interest in that reference comprise 
externally blown flaps attached to the trailing edge of a wing. In the 
present investigation, the computer program of reference 9 was extended 
to handle flaps typical of augmentor wing/flap systems, as shown in 
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figures 2 and 3. Each coroponont is covered with horseshoe vortices on 
its moon surface. The trailing legs of the horsoshoo vortices are wade 
to lie in, and extend to infinity in the plane of each surfoee. However, 
the trailing logs of the horooshoe vortices on the wing (or that surface 
of the ougmontor wing ahead of the forward flap) are hont at the v/ing 
trailing edge so that they lie in the plane of the forward flap as 
illustrated in figure 3. 

The wing/flap geometric parameters and flow conditions to)ion 
account of in tho prosont methods arc suimnarizod below. More detailed 
descriptions are given in roforonco 9. The wing may have 

• finite span 

• breaks in swoop of the loading and trailing edges 

• uniform diliedral over tho span 

• taper sot by loading- and trailing-odgo shapes 

• twist and camber 

• zero wing thic)tnoss only 

Up to 10 flap surfaces can bo accounted for. description of 

the plnnform of the flaps is governed by tho following goomctrlo char- 
acteristics. 

• partial span 

• straight leading- and trailing-odgo shapes over the span 

» taper sot by loading- and trailing-edgo swoops 

• root chord in plane parnllol to vertical or x-z piano, soo 

figure 3 

• tip chord must lie in a vortical plane parallel to the vortical 

plane containing tho root ahord 

• twist and camber 

• zero flap thickness only 

Effects of angle of sideslip and compressibility are not included. Angle 
of attack, flap deflection angle, dihedral angle, twist and camber angles 
are accounted for in the flow tangoncy condition using trigonometric 
functions since some of those angles can be large. 

Vortex layout .- A few horseshoe vortices and the reference coordinate 
system (x,y, z) are shown in figure 3 for a tapered augmentor wing/flap 
system with swept leading- and trailing-edgcs. Mie flaps are set at 
large deflection angles typical of hover and transition flight. An 
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aygmontior diffuaor soction is for.tiod by the foKward and aft flaps. Tho 
chardal pianos of the wing^ center flap# and diffusor flaps are divided 
into trapoKoidal area olemants or panels. On the forward flap, the 
spanwiso distribution of the panels must bo the Mmu as on the wing. A 
horseshoe vortex is placed in each panel such that the bound (or spanwiue) 
log lies along the panel quarter ehord and the trailing legs lie along 
the side edges of the tmnel. a’ho trailing logo extend back to infinity 
in the plane of the panel except for the horseshoe vortices on the wings. 
The legs originating on this surface trail back to the trailing edge 
and are then bent to lie in the plane of tho forward flap. In fact, the 
wing trailing logs will coincide with the trailing legs of tho horseshoe 
vortices on tho forward flap. Figure 3 shows one horsoshoo vortex only 
on the chordal planes of tho wingi forward flap, center flop and aft 
flap. On the forward flap the trailing logs associated with tho wing 
vortices are shown slightly separated from the trailing legs of the 
horsoshoo vortices of the forward flap for tho purpose of clarity. The 
strengths of the horsoshoo vortices, unknown so far, aro obtained from 
tho flow tangoncy condition doacribed next, 

Boundary condition.- The flow tangency boundary condition is applied 
at control points located at tho midspan of tho throe-quartar chord lino 
of each area element or panel. Details are given in roferenco 9, pages 
9 through 14. A brief account only is given here. Some control points 
are indicated on figure 3. Dosignato tho total resultant velocity 
vector at one control point as q and tho normal to tho surfaco in 
question as h, then tho boundary condition states 

q • n « 0 (1) 

for the finite set of control points, with no jot present, velocity q 
includes contributions from all the horseshoe vortices laid out over the 
augmentor wing surfaces and tho free stream component. The velocity 
components induced by the horseshoe vortices arc related to tho unknown 
vortex strengths through tho use of influence functions given in rofer- 
onco 9. The functions associated with horseshoe vortices on the wing 
account for tho effects of the angle of deflection of the trailing legs. 
In the application to augmentor wings, velocity contributions induced by 
the jot singularities are included in equation (1) for tlio purpose of 
accounting for jet interference. A later section concerned with the jet 
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inoaol doieribos tho mothod ugod to gonorate the jet AnducGd volocity 
contribubionSf Not© that tho fuoe stironm componont roproaontg tho 
flight conditions hovor (goto volocity ) , transition or forward flight 
(nonzero) , 

'fho application of tho flow tangency condition/ oquation (1 ) , at 
tho control points distributed over tho chordal pianos of tho ouymontor 
wing surfaces rosulbo in a sot of simultv'inaous equations from which tho 
unknown vortoK strengths are obtained, Once the strengths are known, 
the aerodynamic forces and morionts acting on the augwontor wing compo- 
nents aro calculated using tho method doacribod next. At this stage it 
is also possible to compute, at any field point, the flow field induced 
by the horseshoe vortices representing tho surfaces of tho ougmontor 
wing. 


Forces and moments . ^ Tho aorodynamie forces acting on one olomontal 
area or panel can be determined from tho application of tho Kutta- 
Joukowski law for forces noting on a vortox filament. Tho fundamental 
statomonUn are given in reference 12. This approach has boon used 
successfully in reforancos 9 and 10 and others, Details of tho applica- 
tion of the vortex filament force method are discusr.od in reference 9, 
pages 17 and 10. A sununarizod account only will be given hero. 

The aerodynamic force actincj on a vorto.x filament per unit length 
of filament is expressed as the vector product of the velocity of tlie 
flow q past tho vortex of strength r. 

F « q X r (2) 

Contributions from all horseshoe vortices laid out over the augmentor 
wing components, free stream and jot induced contributions are included 
in the calculation of q. The elemental panel force calculated with 
equation (2) can be used to compute spanwiso load distributions by sum- 
ming over all the panels and overall pitching moments are determined by 
using the panel forces times the appropriate moment arms. Tlie results 
include normal force, axial force and pitching moment coefficients. Tho 
directions of the normal and axial force coefficients, and C^, 
respectively, are shown in figure 3. 
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Jot Kodol Including Jot Wako 


Tho nnalyfcicnl singula;city distribution chosen to ropresent the 
extornal induced o££octs of the jot is doscribod in this section. I'he 
jot, exhausting from the trailing odgo of tho contor flap (figure 3 ) , 
is divided into two distinct regions for modeling. Tho first region i» 
that portion of tho jot from its exhaust notxlo to the exit of tho dif- 
fusor, Secondary jots arc not included in the prosont analysis, Tho 
second region is tho remainder of tho job from the diffuser exit to its 
chosen end point at a finite distance downstream of tho exit. Tho 
following is a description of tho details of this jet model. 


Gonoral approach . ~ Tho objocUivo of tho job model is be represent 
the induced external flow of foots of the jet, ignoring tho details of 
the mixing flow inside the diffusor and tho internal flow of the initial 
region of the jot. Tho basic flow model used as a starting point for tho 
required jet model is that of an actuator disk which can bo used to 
roprusont a jot of fluid with Jjigher velocity and higher total hood than 
the surrounding fluid. Thrs flow model has been used to represent the 
external flow field induced by tho v/ake from a turbo jot or turbofan 
engine (rofo. 9 and 10). In roforonce 9, the development of tho flow 
model for circular and elliptic cross section jots is described, and in 
reference 10, tho extension of this model to a rectangular cross section 
jot is prosontod. Tho latter rectangular jet model is tho most appropri- 
ate model for tho application to augmontor wings. 


The equivalent singularity distribution for an actuator disk is a 
semi-infinite length cylinder with a uniform distribution of vorticity 
on its surface (ref. 13) . Two characteristics of this model are a 
uniform velocity profile inside the job and an increasing mass flow with- 
in the boundary over the initial few diameters of length. This model 
has boon extended to expanding boundaries, various cross sectional shapes, 
and bent centerlines for use as a turbojet wake model. In making those 
modifications, the analytical singularity model of a uniform vorticity 
distribution on a semi-infinito cylinder vvas changed to a finite length 
distribution of vortex rings. This cxiange was necossitatad by numerical 
difficulties in calculating tho induced velocity field associated with 
the singularity distribution for curved jobs, and more details of this 
change are presented in reference 10. 
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The above mentioned singularity distribution v?ill be adapted to the 
current problem of developing a jet model for an augmentor wing. Pown-* 
stream of the augmentor, the jot should behave as a free jet and thus 
can bo handled in tho same manner as the previous jot wake models. The 
upstream region of the jot model, in tho diffusor region, is quite dif- 
ferent from previous jot models. For tho case at hand, tho jot o.xhausts 
from a very high aspect ratio slot nosslc at the trailing edge of tho 
center flap and expands very rapidly to fill the diffusor. In this 
process tho jot entrains air drawn into the diffusor, Tho assumption 
is made that tho jot and entrained air are completely mixed and fill tho 
diffusor exit. The rates of expansion and entrainment of secondary 
fluid by the jet are much larger than typical free jot rato,^ bocause of 
the enhanced mixing over the relatively short length of the diffusor. 
Therefore, a now flow model is required in this region. This is described 
in tho following sections , 

Diffusor region .- This portion of the jet model, between the center 
flap trailing edge exhaust nozzle and the end of the diffusor section, 
is quite different from the jet models described in references 9 and 10. 
The major difference is the high rate of expansion of the jet boundaries 
due to the mixing in the diffusor. In actual practice, the mixing process 
is enhanced by tho use of hyper-mixing nozzles and boundary layer blowing 
on the inner surfaces of the diffusor. The entrainment rate of tho jet 

is therefore much higher than that of a free jet because of tho mass 

augmentation effect of the diffusor. In this context, tho fluid entrained 
by the jet is assumed to be the entire amount of secondary fluid; that 
is, the entrained flow is the difference between the total mass flow at 

tho diffusor exit and the jet mass flow at the primary nozzle. 

The boundaries of the jet in the diffusor region are shown dashed 
in figures 3 (and 6). Since the present model is not concerned with the 
details of the internal flow, the idealized boundaries are prescribed by 
straight line segments as shown. The only requirement is that the jet 
model expand to fill the exit at the end of the diffusor. The initial 
portion of the jet is sized to match the actual cross sectional area of 
the jet nozzle at the center flap. The short length of nonexpanding 
boundaries near the trailing edge of the center flap, shown in figures 
3 (and 6) is included to give the analytical jet model an initial run 
length to build up to the correct jet velocity. This length is typically 
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foux* or five timos bho minimum cross soction dimension, in this case, 
the v/iduh of tho jot. sinco dotails of fcho actual jot in el'iG diffuser 
arc unknown, tho jot boundary is specified to expand linearly to the 
exit of the diffuser after the initial run length, and tho centerline 
of the jet is positioned approximately midway between the forward and 
aft flaps. The spanwiso dimension of UJia jot is assumed constant; 
therefore, the width of tho jet is the only changing dimension. This 
results in a dooroasing eross-soction aspect-ratio along tho length of 
the jet in tho diffusor. A schematic of tho jot model, to bo u.sod later, 
is shown in figure 4 (a) . 

Tho distribution of vorticity, along tho length of the jot 
boundary is determined in the following manner. Tho initial vorticity 
.ntrongth nt the primai'y jet noxxlo is specif iod from tho known value of 
primary nozalo mass flow velocity v.. to bo 


This vorticity strangth per unit length remains constant in the initial 
region to give tho centerline velocity an opportunity to stabilise at the 
correct value. Since the jot model does not attempt to model the dis- 
tribution of velocity inside tho diffusor, the next point at which jet 
mass flow is known is at the diffusor exit. At this station the vor- 
ticity strength por unit length is specified to bo 


Ye “ V. (4) 

where v.\ is the flov; velocity associated with the total mass flow at 
“'h 

tile end of tho diffusor. 

Tho distribution of vorticity between tho initial and exit values 
was examined using two different mothods. The first method is a linear 
variation of strength between the two known values. The local vorticity 
strength is given by the relation 



(5) 


as illustrated in figure 4 (b) . This disteibution was abandoned after 
some preliminary calculations. Tho discontinuity in tho slope of the 
strength curve at each end caused certain difficulties in the calculation 
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of the inducocl volocitios insUlo fcho jot boundary. Tho induced vcloc- 
itioB on tho contoriino could not stabilisic at the appropriate initial 
and Qxit values bocauoe of tho rapid chango in vortox strong tli ooqurring 
at tho boginning and ond of the linoar region. 


I'ho second approach to tho apoeifieation of tho vorticity botwoon 
and Yg uses a ninth”ordor polynomial distribution. It is dofinod as 


9 8 7 6 5 4 

« ar + hfi + of,' + df, 4- of. + f| 

^13 


+ gf,^ + hf,^ + if. 


( 0 ) 


where 


'E 


B 


Ye 


(7) 


with rospect to f. 


We require that tho first tliroo dorivativos of _ 

Yf " Yj .. 

be tore at S *- and tho first si.x derivatives be zero at s Sj|,. in 
this way, the vorticity is concentrated towards tho jet exhaust. 'Jl'his 
is necessary to model tho liigh jet-flow volocity at the exhaust nozslo 
and allows for high ontraimnenb rate.s ovor the length of the diffusor. 
Equation (G) becomes 


Y *" T p o o 7 

^ « 28f,- - 36C (8) 

~ 'e 

which produces smooth conterlino velocity distributions between and 
for any values of Yj. and Yj^* 

conterlino velocity distributions for two thrust levels are shown 
in figure 5. From tho end of tho initial region of the jot to tho end 
of the diffusor region, the velocity is very smooth and wall behaved. 

The predicted volocity in the initial region is typical of the results 
obtained with this jot model. 

One coamient can bo includod at this time regarding tho possibility 
of calGulating the strength of the vorticity distribution which will 
pi*oduco a specified conterlino velocity distribution. 'Wiis ai^proach was 
axaminod with partial success. If the expansion of the jet is moderate , 
for example, if bg « -Ib^. (b is the jet width), then this inverse pro- 
cedure works very well. It does have the disadvantage that velocities 
must bo specified at a large number* of points on the conterlino, but this 
can bo handled by curve fitting or assuming a distribution of volocitios. 
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Whan tho expansion is much larger, numerical difficulties arise which 
prevent a reasonable solution for the vorticity strengths. Those dif- 
ficulties stem from an ill conditioned matrix in which the diagonal 
terms arc not dominant. The results indicate an erratic distribution 
of vorticity which do not produce a realistic external velocity field. 
For this reason, this approach was abandoned in favor of the method 
employing the ninth-order polynomial distribution. 

Downs troam region .- The portion of the jot model in the region 
downstream of the diffusor e.xit is treated as a free jet wake. The 
length of the jot is chosen to be sufficiently long so that any further 
increase in length does not affect appreciably the flow velocities pro- 
dictod by the jot model at the exhaust noxsle and the diffusor exit. 

The vorticity distribution on the wake boundary is hold constant at the 
value specified at the end of the diffusor, which is the roquix‘oment 
for a free jet model (ref. 10) . The remainder of the jot is specified 
in a manner similar to the procedure presented in roforenco 11. 

The wake cross section is rectangular over the entire length of the 
wake, but the boundary is allowed to expand according to available 
empirical information. Holding the span of the wake constant, the 
spreading information contained in roforenco 11 is applied to the wake 
width. The effect of an expanding boundary does not have a largo effect 
on the induced loading on the wing and flap surfaces; therefore, this 
approximate approach for sotting the spreading rate is considered 
adequu'Q for the current investigation. 

The path taken by the jet after it loaves the end of the diffusor 
must be specified with respect to the location of the diffvisor exit. 

This is done using empirical information on the patJi of a rectangular jet 
in a crossflow (ref. I*!) . Small difforoncos in the contorline path do 
not cr'Cato largo difforoncos in the induced loading on the lifting sur- 
faces; thereforo, based on the success of a similar approach in refer- 
ence 11, the empirical data of reference 14 are used for aJ.l calculations 
included in this report. The jet wake boundaries are shown in figure Ca 
tor tho oaso - 4 and in tignro 6b for .. 1.4. 

The final component of the jet wake model in the downstream region 
is the blockage model. It is well known that a jet exhausting into a 
crossflow behaves as if the jet boundary is nearly a solid surface. To 
approximate this effect, tho surface of the wake is represented by a 
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finite munbor of vortox quaclriXaUoral panels with a control point at tho 
panel controid as shown in figure 6. The boundary condition of no flow 
through tho control point on each panel results in a set of simultaneous 
equations. The velocity to bo canceled at the blockage panel control 
points consists of a contribution of the free stream to which are added 
the perturbation velocity con\pononts induced by tho vortex lattice on 
the wing/fiap system with power off. This has tho effect of allowing 
the blockage panels to bo porous to tho fluid entrained by tho distribu" 
tion of vorticity modeling the wake (as well vns the jot) . It also takes 
account, in the first approximation, of the effects caused by the aug- 
mentor wing on tho jot wake. The blockage panel or quadrilateral vortex 
strengths are then determined from the simultaneous equations. Vcloc- 
itios induced by tho blockage panels will be included in tho calculation 
of tho powor-on horseshoe vortex strengths of the wing/flap system at a 
later stage. 


Method of Solution 


The flow models discussed above have boon implemented in computer 
programs which arc arranged to bo used in a sequential manner. Tiiis 
arrangement was found to be convenient in the iteration scheme used to 
arrive at a solution. In a series of stops, tho vortex lattice and jet 
analysis arc applied as follows: 

Stop 1. The vortex lattice analysis is applied to the augmentor 

wing surfaces without jot induced effects in the boundary 
condition. The power-off horseshoe vortex strengths are 
determined. Velocities induced by the horseshoe vortices 
at the control points of the blockage panels are 
computed . 

Stop 2. Tho jet wake centerline and boundaries are located. 

Dlockage panel strengths are calculated using velocities 
induced by the power-off horseshoe vortex lattice witli 
strengths determined in step 1. At this point tho block- 
age portion of tho jet wake model has boon modeled. The 
blockage-induced effects at the control points on the 
wing/flaps surfaces are calculated and will bo used in a 
later stop. 
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stop 3. 


Tho jot model is now applied to tho augmontor jot for 
purposes of calculating the distribution of vorticity 
within tho augmontor and downstream of tho augmontor exit. 
As noted previously, tho vorticity distribution is deter- 
mined by the jot velocity at the exhaust nozzle of tho 
center flap and at the diffusor exit. These are specified 
from tho augmontor performance proscribed initially for 
tho calculation. Tho assumption is made that tlioso 
velocities are produced only by tho distribution of vor- 
ticity reprosonting tho jet, and the vorticity distribu- 
tion is calculated. The jot-induced velocities at tho 
wing/flap control points are calculated. 

Stop 4. The vortex lattice analysis is applied again to tho aug- 
montor wing surfaces. This time, the boundary condition 
includes velocity components induced by tho distribution 
of vorticity and the blockage panels modeling the jet and 
its wake. Tho horseshoe vortex strengths arc recalculated. 
Next,- the flow field at the diffusor exit is computed 
including contributions generated by tho vortex lattice on 
tho wing/flap system, tho distribution of vorticity and tho 
blockage panels modeling the jet, and the component of the 
free stream (zero for tho hover case) . Tho area averaged 
flow velocity is doUorminod and compared with tho specified 
diffusor exit average velocity. 

At this stage, the predicted value is usually higher in magnitude 
than tlic specified value. A lower diffusor exit velocity is selected 
and fed back to the jet analysis, step 3. All other input is kept tho 
same including jet exhaust velocity, centerline location, and spreading 
rates inside and downstream of the diffusor region. For tho same horse- 
shoe vortex lattice and blockage panel layout, step 3 is repeated with 
the adjusted diffusor exit velocity for the jet analysis, A new vortex 
distribution is computed for the jet and the induced velocity components 
at the wing/flap control points ai*e updated. Stop 4 is repeated 
and the average flow velocity at the diffusor exit recomputed. If the 
vipdated value matches tho specified one within a selected error bound, 
the itei‘ation is stopped. The forces and moments calculated by the 
horseshoe vortex analysis now reflect tho effects of the mutual 
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intorforonco botwoon tho jot and tho augmontor wing surfaces and tho 
effect of the free stream. The overall forces and moments are calculated 
as the sum of the contribution from the jot itself (thrust at the exhaust 
noaslo) and tho contribution from forces acting on tho horseshoe vortex 
lattice roprosonting tho wing/flaps system. 

It should bo noted that tho specified velocities at tho exhaust 
nozale (at tho trailing edge of tho center flap) and the diffusor exit 
should bo for tho actual flow conditions at hand. As such, the specified 
diffusor exit velocity deduced from experimental data (or directly 
measured) is already representative of tho presence of tho augmontor wing 
surfaces and especially tho effects of tho actual jot wake. Therefore, 
as long as the iteration scheme described above results in a predicted 
flow velocity at tho diffusor exit that matches the specified value, 
the location and shape of the jot wake is not of primary importance in 
the calculation of tho forces and moments. In other v/ords, under these 
conditions the jot wake has a small effect on loads. Tho region of the 
jet inside the diffusor, through its large entrainment effect, is mainly 
rosponsible for the intorforenco effects of the jot on the aerodynamic 
loads acting on the augmontor wing surfaces. Finally, it is noted that 
this procoduro not only produces tho correct mass flow at tho center flap 
nozzle exit and the diffusor exit, but the correct secondary flow enter- 
ing the diffusor. 

THEOHETICAL RESULTS AND COMPARISONS 

Partially due to tho newness of the augmentor wing concept and 
partially for proprietary reasons, very limited experimental data 
involving augmontor wing configurations are in tho public domain. In 
particular, component loading and pressure distribution data on a VTOL 
fighter- typo configuration such as shown in figure 1 are practically 
unavailable. There is probably more data available for static (zero 
forward velocity) conditions. Even though tho methods described in this 
report can handle the static (or hover) case, it was more important to 
test the methodology with forward flight conditions because of the 
emphasis placed on developing a lifting surface modeling scheme, account- 
ing for jot interference effects, applicable to VTOL fighter- type aug- 
mentor wings in flight. 
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The only data made available to test the augmontor wing analysis do- 
sceibcd heroin involves the wind tunnel transport model of reference 8. 

This model is shown in figure 7. A cross section of its rectangular aug- 
mcntor wing system is indicated in figure 2 for a set of flap dafloction 
angles ropresontative of transition conditions. Inherently/ this model 
does not resemble a VTOL fighter typo, figure 1/ in that the augmentor 
wing system is unswept and thick in section. Compared to the stroamwiso 
flap lengths, the forward or wing part is extremely short. In addition, 
the model employs a fuselage with considerable cross section area atypical 
of a VTOL fighter. Furthermore, figure 2 shows secondary jet nozzlou on 
the forward and aft flap components. In VTOIj fighter configurations, those- 
jets serve to control the boundary layers on the diffusor walls formed by 
the forward and aft flap surfaces. The amount of engine exhaust diverted 
to the secondary nozzles may vary. For the transport model under consid- 
ex-ation, the size of the secondary nozzles in the forward and aft flaps 
indicate that appreciable amounts of exhaust air* could be diverted to thorn. 
References 8 and 15 do not contain information about the division of flow 
between the primary and secondary nozzles. 

In spite of the somewhat unsuitable geometric characteristics of the 
wind tunnel model described above and the uncertainty in the division of 
flow between the primary and secondary jet nozzles, the theoretical methods 
were applied to this configuration for preliminary verification and to 
point out the usefulness of the methods for indicating areas of improve- 
ment in the preliminary design stages of augmentor wings. 

In what follows, the layout of the horseshoe vortices on the surfaces 
of the augmentor wing under consideration is described. Then, based on 
the assumption that the flow out of the secondary nozzles is small com- 
pared to the flow from the primary nozzle, the jet centerline and boundary 
positions are specified and the blockage panel layout is given for two 
power settings. Some of the calculated flow fields induced by the vortex 
lattice alone and jet alone will be shown. They are the results of the 
intermediate steps described in the section entitled "Method of Solution". 
As a result of the iteration, the adjusted diffusor exit velocities input 
to the jet model are given, and final velocity distributions calculated 
by the present method at the exit are shown. Finally, comparisons are 
made between the predicted and measured overall forces and moments. 
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Pnnol Layout on Surfaces 


The horseshoe vortex lattice layout on the surfaces of the augmontor 
wing for use in the comparisons with data obtained with the transport 
model shown in figure 7 will now bo described. As can be soon on that 
figure/ the model includes a fuselage and a tail section. The former 
will only bo partially accounted for by the present method but the latter 
(horizontal tail) is represented by an additional horseshoe vortex 
lattice to be solved simultaneously with the vortex lattice on the sur- 
faces of the augmontor wing. Thus/ as far as the model components are 
concerned/ the augmontor wing system and horizontal tail will be handled 
by the present methods including an approximation for the lift carried 
over on to the fuselage. Plow conditions include zero angle of attack 
and two forward flight speeds. 

Chordal planes .- The wing, forward flap, center flap, and aft flap 
components of the augmentor wing system, sho\-/n in figure 2, are idealized 
to the chordal plane representation shown in figure G. The v/ing/flap 
components arc oxtendod through to the fuselage centerline to account for 
body lift carryover, however, the jet is made to span only over the 
e-xposed wing/flap region and the effects of the jet are felt by tlic sur- 
faces in that region only (i.e., no jet effects are included on the part 
of the augmentor wing surfaces inside the fuselage) . In this way, the 
lift carryover onto the fuselage is accounted for to first order for both 
power-on and power-off conditions. At the present time, the methods can- 
not account for the effects of the fuselage (Beskin upwash) on the aug- 
mentor wing surfaces. Therefore, flow conditions including nonzero angle 
of attack introduce uncertainty with the present method when a fuselage 
is part of the configuration under consideration. 

Symmetry . - It should be noted that because of geometrical and flow 
symmetry about the vertical plane through the fuselage centerline, only 
one half of a given configuration need be covered with a vortex lattice. 
In fact, the transport wind-tunnel model under consideration here is a 
half model mounted on the tunnel wall for the same reasons. 

Number of horseshoe vortices .- For all comparisons discussed later, 
the augmentor wing surfaces on one side of the vertical plane of synm\ctry 
will be treated as follows. In the chordal plane of the wing of the 
rectangular wing/flap system, shown in cross section in figure 6, 11 


horsoshoQ vorticos are laid along the span and three along the chord, On 
the forward flap, the spanwise number is the same but four horseshoe vor- 
tices are laid out on the chord because of the longer length involved. 

On the center flap, the spanwise number is also 11, and throe are placed 
along the chord. The aft flap, with the longest chord length, is covered 
by 11 spanwise and five chordwisc horseshoe vorticos. As a consequence, 
there are 33 horseshoe vorticos on the wing, 44 on the forward flap, 33 
on the center flap and 55 on the aft flap. Note that the trailing logs 
of the Jiorseshoc vortices on the wing are bent at the wing trailing edge 
to lie in the chordal plane of the forward flap. This construction is 
also indicated for the one horseshoe vortex on the wing of the general 
augmentor wing svstem of figure 3, Finally, the horizontal tail is 
covered by an auuitional liorseshoe vortex lattice consisting of five 
vortices along the span and three along the chord. Thus, a total of 180 
horseshoe vortices represent the augmentor wing/flaps and horizontal tail. 
Their strengths are determined from one sot of simultaneous equations as 
described in an earlier section. 

Layout of Jet Boundaries and Blockage Panels, 
and Jot Induced Velocities 

For the case at hand, the jet and jot wake are to be modeled for two 
power settings using the method described above in the section concerned 
with the jot model. The power settings are expressed as the thrust or 
jet momentum coefficient, associated with the jet at the exhaust 
nozzle located at the trailing edge of the center flap. The assumption 
is made that the secondary nozzles have negligible effect. 

Mvj 

C ,= -1- = (9) 

‘^^REF ^®REF 

In the above expression, T is the thrust force produced at the trailing 
edge of the center flap, M is the mass flow rate and, vj^ is the exhaust 
nozzle mass flow velocity referred to earlier in equation (3) . Com- 
parisons with experimental data from reference 8 will be discussed for 
Cjj = 1.5 and 7.31. Note that due to symmetry, the effects of the jet on 
the opposite side of the vertical plane of symmetry must be accounted 
for. At the present time, the mutual interference between the left and 
right jots is not included in the analysis. 
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Layout inside ejiffusor .- Coniiclor tho chordnl plane roprosontation, 
shown in figuiro G, of tho rectangular nugmontor wing of figuros 2 and 7. 

The boundaries of the portion of tho jot inside tho diffusor formed by tho 
forvmrd and aft flaps is indicated by tho dashed lino in figure 6. It is 
based on tho procedure described earlier concerned with tho diffusor region. 
The "oKit" location is defined by the location along tho straight jot 
centerline v/here tho jet boundaries roach their maximum width. This posi~ 
tion is also marked on figuros 4(a), 5 and G. Referring to figures G (a) 
and 6(b), the coordinates of the contorlino and the width b of tho jot in- 
side the diffusor are the same for both power settings and are specified 
below. For tho jet on the loft side of the piano of symmetry, tho geo- 
metrical characteristics are given in tho following table. 
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cm (in) 

cm (in) 

cm (in) 

-20.27 (-7.98) 

-21. G9 (-8.5) 

-29.49 (-11.01) 

-31.01 (-12.21) 

-32.51 (-12. B) 

1 1 
1 

1 ^ O 

68 (“38. 

r ' 

025) 

' 

-3.71 (-1.40) 

-1.22 (-0.48) 

13,07 (5.38) 

16.51 (0.50) 

19.35 (7.02) 

0.30 (0.12) 

0.30 (0.12) 

10.29 (4.05) 

10.29 (4.05) 

10.29 (4.05) 


As can bo seen from the above table and figure 6, the jot boundaries do not 
expand over the initial run length. This is for tho purpose of building up 
the correct primary jet velocity Vj^. Thon, the jet boundary is mode to 
expand linearly to the "exit" of the diffusor. This "exit" location was 
chosen to coincide approximately with tho taps of a total pressuro rako 
used in tho tests described in reference 8. Also note that in this layout 
the forward jet boundary meets the trailing edge of the forward flap. 

In the spanwise direction, the jet and jot wake are positioned from 
tlie side of the body out to the midspan of the chordv/ise row of horseshoe 
vortices at tho tip of the augmettar wing. Xjoenting tho outboard side 
face of the jet (and its wF^ke) inboard from tho side edges in this way 
avoids undue influence at the blockage panels control points induced by 
the strong vorticity along the side edges of the wing and flaps. For the 
case at handr tho spanwise dimension of the jet is 142.37 cm (56.06 in.). 


Input veloclfciog for jot itiodol .- Onco fchss goomstaeical layout of th® 
^Qt Insid® til® diffusor region is dafinod, tho primary ;jQt volocity and 
the volocity of th® fully mixed (augmonted) jot at tho diffusor oxit muat 
ho apocifiod. In addition, the hloehago panels on th® jot wako arc 
exposed to a certain flow field ns will be discussed Inter in this sec- 
tion. For tho two thrust coefficients, tho following table contains tho 
required velocities divided by the free stream. They wore calculated 
using the ideal gas relationship and certain assumptions about pressures 
and temperatures at tho jet nosalo and diffuser "oxit" as oxplainod below. 
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1.50 

7.31 

1.20 

1.48 

3.9 

4.7 

25390 (833) 
32240 (1058) 

4298 (141) 
7407 (243) 

3048 (100) 
1905 (62.5) 

8.33 

16.94 

1,41 

3,89 


At tho primary nozzle exhaust the static pressure was assumed to bo 
1013 . 06.5 newtons/m’" ( 14.7 psia) . Tho total temperature in tho primary 
nozzle flow was assumed to be 15.56°C (60®F) . X^’rom reference 8, onclosuro 
18, the total-to-otatic pressure ratios for tho ijrimary jot are specified 
as 1.5 and 2.0 for C, « 1.5 and 7.31, respectively. The free stream 
dynamic headG associated with these two thrust coefficients are 
569.77 newtons/ro^ ( 11.9 Ibs/ft^) and 221,21 nowtons/m^ (4.62 Ibs/ft^') , 
rcspoctivoly . Density at the diffuser exit was taken equal to free 
stream density. The diffusor oxit velocities, Vj^ (shown above), require 
knowledge of the secondary or entrained mass flow rate and primary 

nozzle flow rate m . X^oference 15, i^agc 687, contains a curve relating 

•!» 

the mass flow ratio, m_/ra^, to the thrust augmentation ratio. Tho thrust 
augmentation ratio, d, is given in roferonco 8, enclosure 20, for tho two 
thrust coefficients listed above. 


The above mentioned mass flow rates are for the actual flow conditions 
at hand and include the effects of jot flow issuing from the primary and 
secondary nozzles. The average diffusor exit velocities arc determined 
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from those data. Thus , even though the analysis assuwos that all the 
exhaust gas is diverted to the primary notzlc, the uses of the above exit 
velocities guarantees that the moaeurod amount of entrainment is 
accounted for in the present method. 

The values given above for the pet velocity at the primary nettle 
exit and the diffuser exit will bo used in the determination of the dis- 
tribution of v-.rticity in accordance with oyuation (8) . As noted earlier 
in the section describing the stepwise method of solution# the diffusor 
exit velocity will require adjustment during the succossiva jet calcula- 
tions in the iteration process. Tho values given in tho above table for 
the two povmr settings will bo used ns the Initial choice. 

Wake contorline locations and nproading rates .- For the case at hand# 
the empirical procedure described in an earlier section under the heading 
"Downstream region" is applied. Tlio jot contarline location is obtained 
for the two diffusor exit velocity ratios shown in tho previous table 
from equation (1) of reference 14# which is based on experimental data. 

The aspect ratio of the jet in a crossflow associated with the cited 
reference is four# whereas the aspect ratio of the augmentor jot is 
approximately seven. Therefore# tho centerline paths obtained for tho 
two power settings using the data of reference 14 are only vipproximate. 
However# in accordance with the arguments given at the end of the sootion 
entitled "Method of Solution"# the coordinates of the jet aft of the dif- 
fusor section need not bo very aecurate as long as tho diffusor exit 
velocity used in the prediction scheme is directly measured or deduced 
from data for tho configuration and flow conditions at hand. 

Tho cross section of the jet wake is rectangular. Proceeding down- 
stream# tho Width of the wake expands but the dimension in the spanwise 
direction is held constant. This latter assumption should be suffi- 
ciently adequate for the large aspect ratio (AR « 7) of the augmentor jet 
under consideration. Reference 11, figure 9# contains information from 
which tho expansion rates arc obtained for tho two diffusor exit veloc- 
ity ratios indicated in tho previous table for the two power settings. 

Note that tho remarks made above with regard to tho accuracy apply to 
the expansion x'ates as well. The following tables contain the cooi‘di- 
nates of tho contcrlino and tho width of tho jot wake for the two power 
settings. Quantity b is tho width of tho jot v/akc. Figure 6 shows the 
coordinate system and the jet wake on the loft side of the plane of 


symmotcy in jaccojrdancc with tho coordinates given in the tables below 
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(23,32) 

17.96 

(7.071) 
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64.54 

(25.41) 

18.04 
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Cji »» 7, 31 
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12.14 
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-45.06 
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(15,63) 
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(5.549) 
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(18.94) 
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80,75 

(31.79) 

22.94 

(9.032) 


For both power settings, the spanv/isc c-xtent of the jot and its wake 
is 142.37 cm (56.06 in.). 
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Blockage panels layout ,- The length of the jet wake (measured along 
the centerline) is taken as approximately three times the chord length of 
the forward flap, see figure 6. Calculations presented in reference 9 in 
connection with externally blown flaps indicate that greater jet lengths 
produce very small changes in aerodynamic loadings at the expense of 
additional computer time. The blockage of the jet boundaries is modeled 
by quadrilateral vortex panels described as follows. 

On the front and rear faces of the jet wake, ten blockage panels are 
laid out in the spanwise direction and seven panels in the lengthwise 
direction. The sides of the jot wake are covered with four blockage 
panels in the direction normal to the centerline and seven panels are laid 
out in the lengthwise direction. Thus, the boundaries of the jet wake on 
the left hand side of the plane of symmetry are covered by 196 panels. To 
preserve symmetry, the same jet wake and its singularity layout is posi- 
tioned on the right hand side of the symmetry plane. 

The velocity field to bo counteracted by the jet wake blockage panels 
is shown in figure 8 in a vertical plane for part of the jet length, in 
accordance with step 1 of the section headed "Method of Solution", the 
flow velocities are generated by the vortex lattice laid out on the sur- 
faces of the augmentor wing and include the free stream velocity. The 
flow vectors indicated in figure 8 arc calculated by the vortex lattice 
described earlier in connection v/ith the horseshoe vortex paneling layout 
(180 horseshoe vortices) on the lifting surfaces. Flow conditions are 
zero degrees angle of attack and zero power setting. Note that the flow 
velocities at the blockage panel control points are directed downwards. 
Since the vortex lattice and blockage panel strengths are solved for 
separately, the flov/ field soon by the jet wake should be generated by 
the horseshoe vortices as influenced by the jet. As a first approxima- 
tion, the power-off horseshoe vortex strengths are used and in the 
succeeding steps of the iterative approach the blockage panel strengths 
are kept constant. This constraint can be relaxed and the flow field 
impressed on the jet wake recalculated once the jet singularities are 
known. At this time, for the sake of economy the former approach is 
adopted. 

Velocities induced by jet model .- In order to provide some under- 
standing of the entrainment properties of the jet, it is instructive 
to discuss briefly the flow field induced by the jet model at points 


on the components of the augmontor wing. For the jet contoriine an<i hound~ 
ary locations given in the above tables^ the induced Flow field in a piano 
parallel to the .\‘“2 plane for the higher (c^ * 7 . 31) power setting is shown 
in figure 9. The direetions of tho flov; vectors at the control points on 
the forward and aft flaps and tho diffusor entrance indicate strong 
inflow or entrainment towards the initial or narrow part of tho jot near 
the trailing edge of tho center flap. In addition, tho flow inside tho 
jet at the diffusor exit is fairly uniform and tho average velocity ratio 
(vjjj/V^) at the diffusor oxit is about 3.8. Note that at this stage tho 
offocts of the horseshoe vortices are not included in the diffusor e.xit 
velocities. 

Predicted velocity Distributions at Diffusor E-xlt 

I’he augmontor wing of tho transport model shown in figure 7 (and 
with more detail in figure 2) is modeled by the paneling layouts and jot 
model specified above. The stepwise procedure described in an earlier 
section entitled "Method of Solution" is then applied. At tlio end of 
the first pass through that procedure, tho flow velocities across tho 
diffusor exit are computed as the sum of tho voi'tcx lattice and jot- 
singularitios-inducod velocities added to tho free stream. For both 
power settings, tho predicted velocities at the exit wore higher than 
tlio specified average flow velocity. In accordance with stop 4 of tho 
solution procedure, tho diffusor exit velocity, serving I'S one of the 
velocity inputs to tho jet model, is reduced and tho jot model rerun 
(stop 3). Tho adjusted jot singularity strongtlis and induced velocities, 
such as shown in figure 9, are then used in step 4. At tho end of stop 
4, the diffusor* exit velocities are calculated again and compared \\?ith 
the specified average flow velocity, \4hon tho predicted average diffusor 
exit velocity matclies (to within a certain margin) the specified one, 
the solution is considered to bo converged. 

For tho Cjj - 1.5 and 7.31 power settings, the following table con- 
tains the final input velocities for the jet model and the average of 
the predicted and specified diffusor exit velocities. Tho last guan- 
titios have been discussed earlier in the section entitled "Layout of 
Jot Boundaries and Blockage Panels, and Jot induced Velocities". 
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Dp 

The calculated flow velocity at the diffusor exit, , is obtained 
from the aroa-avoraged distribution of flow volocitio's shown in figure 
10(a) for » 1.5 and in figure 10(b) for » 7.31, respectively. Note 
that these velocities are gonoratod by the vortex lattice, the jet vor- 
ticity and jet walso bloc)<age panels and include a component of the free 
stream. The velocity vector associated with the free stream is also 
shown. The distribution of flow velocities across the diffusor exit is 
seen to be fairly uniform in accordance with tl30 rassumption made to that 
effect in the description and specification of the jot model. Roforonco 
0 contains some total pressure data obtained with a ra)<e located at the 
diffusor exit for one spanwiso station. The distributions of flow veloc- 
ities deduced from that data are not uniform and in fact show peaks near 
the diffusor walls formed by the inner surfacos of the forv;ard and aft 
flaps. In the tests of reforonco 8, additional jets i.ssue from the 
leading edges and arc directed along the inner walls of tho diffusor for 
boundary layer control. It is possible, therefore, that tho rake measure- 
ments are influenced by these secondary jets which arc not accounted for 
in the present method. 

The agreement between the calculated diffusor exit flow velocity 
and tho specified value shown in the table above is quite good. At 
this (final) stage, tlio forces and moments computed by tho vortex 
lattice method should bo reprosentativo of the augmentor wing loading 
including tl3c effects of the jet. 
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Force and Moment Comparisons 

Lift, drag and moment coefficients are shown as a function of 
thrust or jot momentum coefficient in figures li(a), 11(b), and 11(c), 
respectively, for zero angle of attack. The open symbols are values 
measured on the complete configuration shown in figure 7 for the set of 
flap deflection angles specified. The closed symbols represent calcu- 
lated values for « 0.0, 1.5 and 7.31 obtained with the paneling lay- 
out and jet specification described above for the same flap deflection 
angles. Only tlie augmentor wing and horizontal tail surfaces are 
accounted for in the present mot)iod; thoroforo, any effects from the 
fuselage on the aerodynamic loads are not included in the calculated 
values. As far as the drag is concerned, the theory calculates the 
induced (due to lift) drag contribution only. 

At the zero and lower (C^ = 0.0, 1.5) power settings, the agreement 
between moasuromont and theory is good. At the highest moment coeffi- 
cient (C^ « 7.31) the lift and drag coefficients arc overpredicted and 
the pitching moment is underprodicted. Overall, the predictions show 
similar trends as the experiment foe increasing thrust coefficient. 

Reference 8 contains a few chordwiso pressure distributions for one 
spanwise location for the selected flap settings. Especially at tlie 
higher power setting, flow separation is indicated on the center and aft 
flap. The experimental pressure distribution along the upper surface of 
the wing and forward flap components (see figure 2) indicate strong 
suction pressures aft of the secondary nozzle. Similar behavior is 
indicated at the loading edge of the aft flap on the surface ahead of 
the secondary nozzle. These observations and the nonuniform measured 
total pressure across the diffusor exit discussed earlier seem to indi- 
cate strong blowing out of the secondary nozzles. This behavior makes 
itself felt more strongly at the higher power setting since the total 
pressure distribution at the diffusor exit is shown to be more uniform 
for the Cjj - 1.5 case. Thus, the partial separation and the strong 
effects of the secondary nozzles at the higher power sotting (C^j = 7.31) 
may account for the discrepancy bet\>;een theory and experiment for that 
condition. 

When interpreting the aerodynamic forces acting on the individual 
flap surfaces of the augmentor wing system, it is often helpful to 
analyze the flow field in which those components are imn\ersed. The 
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present method is capable of determining those flow fields. If, for a 
given wing/flap segment, the oncoming flow field is at a large angle to 
it, partial flow separation over that surface would be expected to occur 
and subsequent loss of lift may result. 

As an o-xamplo of a flow field around the augmentor wing at the same 
flap deflection angles, consider figure 12. Flow vectors are shown for 
the case of aero power. It is seen that the wing is subjected to an 
average upwash of about 15“ . The center flap is immersed in a flow field 
at about 90“ to it, and is likely to suffer from flow separation and its 
consoqvioncos on "lift", it is interesting to note that the aft flap is 
immersed in flow largely aligned with it thus reducing its lift effoo- 
tivenoss. Upstream of the forward flap, there appears to be a low veloc- 
ity recirculation region. Note that at ns much as two chord lengths 
upstr-oam, stx-ong interforcnce effects are predicted by the present method. 

CONCLUDING REMARKS 

A method has been developed for determining the external aero- 
dynamics of VTOL fightor-type augmentor wings in liover or transition 
fliglit for prescribed ejector jot characteristics. Specifically, the 
flow velocity at the primary nozzle and an average velocity at tlie dif- 
fuser exit must be proscribed for the flight condition at hand. The 
augmentor wing may have swoop and taper. The method is based on poten- 
tial flow theory and attached flow is assumed. 

Basically, liorseshoo vertex lattice thoory is used to roprosont the 
solid surfaces of the augmentor wing components. Effects of the jot are 
included in the flow tangency condition. The jet is modeled by a dis- 
tribution of vorticity on its boundary. In addition, blockage panels 
are placed on the jot boundaries downstream of the diffusor exit. Inside 
the diffusor, the jot model is made to expand from the px'xmary nozzle 
width to the full width of the diffusoar exit as seen in sideview. The 
jet model accounts for entraiaamont of secondary air. For given primary 
nozzle and diffusor exit jet velocities, the aerodynamic loadings acting 
on the augmentor wing arc obtained as the result of an iteration scheme 
which produces a calculated average diffusor exit flow velocity to match 
the specified one in magnitude. The iteration scheme consists of a 
serial application of the vortex lattice and jet model methods. 
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Exporimontal data suitable for tasting the developed thoorotical 
methods arc scarca and/or very difficult to acquire. Some comparisons 
with e.xperimontal data taken with a v/STOL transport model equipped with 
an unswopt and untaporod augmentor wing arc shown. This configuration 
is not representative of a VTOL fighter in that the stroamwise sections 
of the augmentor wing components are very thick and it appears that the 
division of exhaust gas flow between the primary nozzle and secondary 
nozzles is not representative. At the present time, the method can only 
accoimnodatc a primary jet issuing from the trailing edge of the center 
flap located above the diffusor. Consequently, the assumption was made 
that all thQ engine exhaust is fed to the primary nozzle. The actual 
transport model was idealized to chordal plane representations of the 
augmentor wing components and horizontal tail. Effects of the fuselage 
are accounted in part only. With those simplifications, the present 
method calculates lift, drag and pitching moments which agreed well V'^ith 
the experimental data at low power settings. For the higher power 
setting, the lift and drag are overestimated and the pitching moment is 
underestimated; however, the predictions show the same trends as experi- 
ment for increasing thrust. 

With all the engine exhaust assumed to emanate from the primary 
nozzles, the total entrainment by the single jot is made to be the same 
as in the case when the power is divided into primary and secondary 
nozzles. As a result, the calculated overall forces and moments should 
be representative but detailed information such as component loadings 
require more knowledge about the flov/ division. 

Notwithstanding the above mentioned limitations, the usefulness of 
the present method includes the capability of indicating potential prob- 
lem areas during preliminary design of an augmentor wing system. The 
developed methods can be used to map the flow field in the vicinity of 
the augmentor wing for determining intorference on other components of 
the aircraft. In addition, the flow field impressed on each of the flap 
surfaces of the augmentor wing system can bo analyzed so that secondary 
nozzles for boundary control can be positioned and sized on the basis of 
that knowledge. Furthermore, the present method can be applied to an 
existing configuration for which component loads have been measured. By 
comparing the component load prediction with measurement, components 
suffering from flow separation and stall can be identified and remedies 
effected. 
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on the basis of tho work performed so far, tho following roeoimen- 
dntions arc offered. 

1. To validate the present method further, additional comparisons 
should be made v^ith available data for different flap settings 
and nonsoro angle of attack. Component load comparisons should 
be made. The majority of tho engine exhaust should bo applied 
to the primary nosslos for comparison purposes. 

2. Detailed data should be taken with a XPV~12A type augmentor 
wing and made availablo for testing tlie developed method. 

3. based on the outeomo of 1 and 2, the single augmentor wing 
system model can be improved by studying effects of different 
jet boundary layouts and jot vorticity distributions (affecting 
ontx'ainraent) inside tho diffusor. Also, tho effects of tho jot 
wake on overall and component loads should be determined. 

4 . It is possible to circumvent the specification of experimentally 
deduced jot voloeitios at tho primary nozslo and tho diffusor 
exit. Detailed internal flow analyses have boon dovolopod else- 
where (ref. 6) capable of generating the required quantities 
for the jot model of the present method. 

5. The present method can be extended to account for more than one 
jet, i.o. to handle secondary jets and aecount for tho associ- 
ated Coanda effects on the aerodynamic loads acting on the aug- 
mentor wing components. Tho applicable technology has been 
dovolopod in connection with USB (Upper Surface Blowing) work 
described in roforence 10. 

6. The present method can be extended to account for augmentor 
canard/augmentor wing systems attached to a fuselage, accounting 
for mutual interference between the canard, wing, and fuselage. 

Nielsen Enginooi*ing & Research, Inc. 

Mountain View, California 94043 
December 1978 
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Figure 2.- Chor<awise 3ection of transport model 
augmentor wing of reference 8. 
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Figure 3,- Chordal planes of an auginentor wing with 
one horseshoe vortex shown on each surface. 
Idealised jet centerline and houndarios. 
Reference coordinate system. 
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Figure 6.- Conclvided. 
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Figure 8.- Volocities induced by vortex lattice on 
augmentor wing surfaces at control points of 
blockage panels on upper part of jet wake at 
y *= -103.71 cm (-40.83 in.) 
spanwise location. 
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Figure 9.- Velocities induced 
points on surfaces of tlie 
y = -77.85 cm (-30.65 


by jot and jet wake at control 
transport augmentor wing at 
in. ) spanwise location. 
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Figure 10.— Concluded. 
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(c) Pitching moment. 


Figure 11.- Concluded. 






Figure 12,- Calculated flow power-off field in the vicinity of the 
augiEeator wing at spanwise location y = —74.93 tm {-29.5 in.). 










